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ABSTRACT 


The  problem  of  interfacing  older  military  terminal 
equipments,  which  work  exclusively  from  an  internal  timing 
source,  into  a  communication  system  that  requires  all  terminals 
to  synchronise  to  the  network  clock  is  addressed  in  this  paper. 
These  data  communication  systems  would  typically  operate  with  a 
fixed  symbol  transmission  rate  and  employ  end-to-end 
convolutional  coding  for  channel  errors.  A  scheme  is  proposed 
which  would  allow  a  clock  discrepancy  of  up  to  +4%  between  the 
terminal  clock  and  the  network  clock.  The  method  proposed  is  to 
add  control  bits  to  the  information  stream  between  terminals. 
The  control  bits  would  allow  the  effective  information  rate  to 
either  increase  or  decrease  by  ±4%.  The  overall  symbol  rate  of 
the  communication  link  is  held  constant  by  selecting  a  higher 
rate  convolutional  code  from  the  set  of  rate  compatible 
punctured  codes  (RCPCs).  This  paper  describes  the  methodology 
used  to  select  an  optimal  RCPC  and  presents  some  examples.  The 
probability  of  error  for  a  selected  code  is  compared  with  that 
of  a  convolutional  code  for  an  application  typical  of  a 
satellite  communication  system.  The  conclusion  is  that  the 
increase  in  signal  to  noise  required  to  sustain  the  same 
probability  of  error  is  minimal. 


f 


iii 


RESUME 


Le  probleme  d' interface  entre  de  vieux  terminaux 
militaires,  fonctionnant  exclusivement  a  l'aide  d'une  horloge 
interne,  et  un  systeme  de  communication,  requerant  la 
synchronisation  de  tout  terminal  avec  1 'horloge  du  reseau,  est 
presente  dans  cet  article.  Ces  systemes  de  communication  de 
donnees  operent  typiquement  a  un  taux  fixe  de  transmission  de 
symboles  et  emploient  un  codage  convolutionel  a  chaque 
extremite  pour  les  erreurs  de  transmission.  Un  arrangement 
permettant  un  ecart  maximal  de  ±4%  entre  1 'horloge  du  reseau  et 
celle  du  terminal  est  propose.  La  methode  suggeree  consiste  a 
ajouter  des  bits  de  controle  au  debut  d ' inf ormation  entre  les 
terminaux.  Ces  bits  de  controle  permettraient  une  augumentation 
ou  une  diminution  de  4%  du  taux  effectif  de  transmission 
d' information.  Le  taux  global  de  transmission  de  symbole  de  la 
liaison  en  communication  est  maintenu  constant  en  selectionnant 
un  code  convolutionnel  a  taux  plus  eleve  a  partir  d'un  ensemble 
de  codes  perfores  ayant  un  taux  compatible.  Cet  article  decrit 
la  methodologie  utilisee  pour  selectionner  un  code  perfore 
optimum  et  presente  quelques  exemples.  La  probability  d'erreur 
d'un  code  choisi  est  comparee  a  celle  d'un  code  convolutionnel 
pour  une  application  typique  d'un  systeme  de  communication  par 
satellite.  L' article  conclut  que  1 ' augmentation  du  rapport 
signal-bruit,  requise  pour  obtenir  la  meme  probability 
d'erreur,  est  minimale. 
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EXECUTIVE  SUMMARY 


Many  modern  Military  Satellite  Communication 
(MILSATCOM)  Systems  assume  that  user  equipment  will  be  able  to 
synchronise  to  a  network  clock.  However,  this  is  often  not 
feasible  with  older  military  equipments  which  have  no  input  for 
an  external  clock.  To  solve  the  possible  rate  discrepancy  the 
system  interface  to  the  equipment  must  either  speed  up  or  slow 
down  the  overall  end-to-end  data  rate.  This  paper  proposes  a 
solution  to  this  problem  for  data  communication  systems  with 
fixed  symbol  transmission  rates  and  convolutional  coding  for 
channel  errors. 

The  usual  forward  error  correction  (FEC)  code 
specified  in  the  MILSATCOM  systems  is  a  rate  1/2  convolutional 
code.  If  a  slightly  higher  rate  code  could  be  found  (i.e. 
13/24),  extra  overhead  bits  could  be  added  to  the  data  stream 
without  impacting  the  coded  data  rate.  The  overhead  bits  would 
be  used  to  indicate  that  either  more  than  normal  data  is  being 
sent  or  less  then  normal  data  is  being  sent,  and  to  carry  the 
additional  data. 

This  paper  uses  a  rate  1/2  convolutional  code  and 
then  punctures  it  by  deleting  parity  check  symbols  to  obtain  a 
rate  13/24,  rate  compatible  punctured  code  (RCPC).  This  code 
permits  the  inclusion  of  2  overhead  bits  for  every  24  data  bits 
transmitted.  The  puncturing  method  chosen  is  optimised  for 
minimum  impact  on  error  rate.  Using  error  bounding  and 
simulation,  the  RCPC  code  is  compared  to  the  convolutional  code 
for  a  typical  MILSATCOM  application.  The  result  is  that  to 
support  a  probability  of  error  of  10~  ,  the  signal  to  noise 
ratio  required  for  the  RCPC  is  0.3  dB  above  the  convolutional 
code.  The  scheme  using  the  RCPC  code  will  support  a  clock 
discrepancy  of  +4%  between  the  data  source  and  the  network. 
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1 .  INTRODUCTION 


Many  existing  military  communications  equipment  do 
not  support  external  clocking,  but  work  exclusively  from  an 
internal  timing  source.  If  the  equipment  is  connected  to  a 
network  and  its  internal  clock  is  not  exactly  matched  to  the 
network  clock,  then  the  interface  will  pass  too  many  or  too  few 
data  bits  per  second.  The  end  result  will  be  a  communication 
link  that  regularly  loses  information  bits  and  synchronisation. 
Any  data  communication  system  that  must  support  this  type  of 
interface  must  rationalise  the  inconsistency  in  clock  rates. 
Unfortunately  many  networks  assume  that  all  connected 
equipments  accept  an  external  clocking  source  or  have  an 
internal  clock  that  is  matched  exactly  to  the  network  clock. 
This  document  proposes  the  use  of  a  rate  compatible  punctured 
convolutional  (RCPC)  code  in  the  data  link  standard  of  the 
communication  systems  to  support  equipment  that  works 
exclusively  from  an  internal  clocking  source. 

The  data  interface  problem  is  described  in  more 
detail  in  Section  2  along  with  a  general  concept  for  correcting 
the  problem  using  a  variable  rate  code.  Section  3  defines  what 
a  RCPC  code  is  and  the  specific  code  parameters  used  in  this 
document.  Section  4  outlines  how  the  optimum  puncturing  table 
is  derived.  The  simulation  of  a  communications  system  model 
used  to  verify  the  theory  is  presented  in  Section  5.  Section  6 
gives  the  results.  Enhancements  to  the  scheme  to  counteract 
jamming  are  given  in  Section  7  and  other  uses  of  RCPC  are 
proposed  in  Section  8.  Section  9  presents  the  conclusions  of 
the  paper. 
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2.  STATEMENT  OP  PROBLEM 


Many  data  communication  systems  operate  with  a  fixed 
symbol  transmission  rate  and  end-to-end  forward  error 
correction  (FEC)  for  channel  errors.  These  systems  typically 
assume  that  the  user  equipment  would  be  able  to  synchronise  to 
the  network  clock  in  the  interface  equipment.  However,  this  is 
not  possible  with  many  military  equipments  which  have  no  input 
for  an  external  clock.  To  solve  this  problem  the  interface 
equipment  must  either  speed  up  or  slow  down  the  effective 
end-to-end  information  rate.  The  use  of  an  elastic  buffer  will 
only  provide  a  temporary  solution  to  the  problem  and  eventually 
an  overflow  or  underflow  will  occur  (depending  on  whether  the 
input  rate  is  too  fast  or  too  slow) . 

The  FEC  code  typically  used  in  MILSATCOM  equipment  is 
a  rate  1/2  convolutional  code.  If  a  slightly  higher  rate  code 
could  be  found  (e.g.  13/24),  extra  overhead  bits  could  be  added 
to  the  data  stream  without  changing  the  coded  data  rate.  The 
overhead  bits  would  carry  any  data  overflow  and  would  be  used 
to  indicate  more  information  is  being  sent  or  less  information 
is  being  sent  than  normal. 

For  consideration  in  this  paper,  data  is  assumed  to 
be  normally  blocked  into  groups  of  24  bits  and  then  encoded 
into  48  bits  with  a  rate  1/2  encoder.  If  by  using  the  higher 
code  rate  two  extra  bits  were  added,  the  network  would  then  be 
able  to  support  a  variation  of  +1  bit  in  24,  or  a  user  clock 
rate  discrepancy  of  t*%.  In  order  for  the  overall  coded  data 
rate  to  remain  constant,  a  rate  26/48  (or  13/24)  encoder  will 
be  required.  Table  1  describes  one  method  of  using  the  two 
extra  bits  to  identify  the  number  of  data  bits  in  the  block  and 
carry  the  additional  data  bit  when  it  occurs.  If  a  constant  24 
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bits  of  data  were  being  transmitted,  the  logic  would  alternate 
between  sending  23  bits  and  25  bits. 


BIT  LOCATION 

RESULT 

1 

2  3  ...  26 

0 

unused  |  data 

23  data  bits 

1 

data 

25  data  bits 

Table  1.  Definition  of  overhead  bits  in  data  link  standard. 


3.  DEFINITION  OF  RCPC  CODES 


The  use  of  convolutional  codes  usually  means  that  the 
coding  rate  is  fixed  to  simple  ratios  (1/2,  1/3,  2/3). 
Non-simple  ratios  result  in  complicated  encoders  and  decoders 
which  are  difficult  to  implement;  however  another  option  exists 
in  the  form  of  punctured  convolutional  codes.  The  concept  of 
punctured  codes  was  first  introduced  by  Cain,  Clark,  and  Geist 

[1]  to  obtain  a  simpler  Viterbi  decoding  scheme  for  rate  2/3 
and  rate  3/4  codes.  This  concept  was  then  expanded  by  Hagenauer 

[2]  to  make  the  coding  scheme  fit  the  transmission  rate.  The 
general  concept  of  a  RCPC  code  is  to  puncture  a  low  rate  1/N 
code  with  period  P  to  obtain  a  family  of  codes  with  rate 

P/(  P+U ,  where  C  can  be  varied  from  1  to  (F-l)P.  From  this 
family  of  codes,  a  code  is  derived  with  the  appropriate  code 
rate.  To  obtain  the  rate  13/24  code  required  for  the  problem 
defined  in  this  paper  the  following  parameters  were  used 

N  -  2,  initial  code  rate  -  1/2 
P  -  13,  C  can  be  varied  from  1  to  13 
C  -  11 
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A  punctured  code  can  be  implemented  with  a  scheme 
similar  to  that  depicted  in  Fig.  1.  The  generator  matrix  would 
be  the  same  as  for  any  convolutional  encoder  and  would  operate 
as  per  normal.  The  puncturing  would  be  accomplished  through  the 
use  of  a  puncturing  table.  The  output  from  the  encoder  would  be 
compared  with  the  value  in  the  puncturing  table.  A  zero  in  the 
table  means  that  the  code  symbol  is  not  transmitted,  a  one  in 
the  table  means  that  the  code  is  transmitted.  In  the  example, 
the  4th  bit  of  the  upper  branch  is  not  transmitted  and  the  2nd 
bit  of  the  lower  branch  is  not  transmitted. 

The  decoding  of  the  coded  information  is  relatively 
straightforward.  The  critical  factor  is  that  the  encoder  and 
decoder  must  be  in  perfect  bit  synchronisation.  The  decoder 
would  use  a  standard  Viterbi  algorithm  set  to  decode  the  data 
coded  with  the  corresponding  encoding  matrix.  In  addition,  the 
decoder  would  have  the  same  puncturing  table  as  the 
transmitter.  A  one  in  the  table  indicates  that  the  encoded  bit 
would  have  been  transmitted  and  the  decoder  can  calculate  the 
maximum  likelihood  matrix  as  normal.  A  zero  in  the  the  table 
indicates  that  the  encoded  bit  would  have  been  been  deleted  and 
the  decoder  would  not  be  able  to  use  this  bit  in  the 
calculation  of  the  maximum  likelihood  matrix. 


4.  DETERMINATION  OF  THE  OPTIMAL  PUNCTURING  TABLE 

As  stated  in  [2],  there  is  no  known  constructive 
method  for  determining  the  optimum  combination  of  generator 
matrix  and  puncturing  table  for  a  RCPC  code  family.  Instead  the 
usual  method  of  determination  is  through  a  computer  search,  as 
was  done  in  this  report. 
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Fig.  1.  Simplified  diagram  of  encoder  and  code  puncturer. 


The  following  lists  the  basic  steps  used  to  determine 
the  optimum  puncturing  table. 


1.  Select  a  suitable  generator  matrix,  in  this  report  a 
rate  1/2  code,  constraint  length  7,  with  maximum  free 
distance  was  chosen. 

2.  For  the  value  of  C  selected  ( £  -  11  in  this  report) 
construct  a  puncturing  table. 

3.  Calculate  the  bound  on  the  expected  number  of  bit 
errors  for  the  generator  matrix  and  puncturing  table 
selected.  The  method  to  calculate  the  upper  bound  was 
first  proposed  by  Viterbi  and  Omura  [3]  and  is  given 
in  Appendix  A. 

4.  Repeat  steps  2  and  3  for  all  possible  constructions 
of  the  puncturing  table. 

5.  Chose  the  puncturing  table  construct  which  yields  the 
best  bit  error  rate  performance  (lowest  upper  bound). 

Using  the  method  described  above  the  following 
parameters  where  determined  for  the  RCPC  code. 


binary 

octal 

Generator  matrix 

1111001 

171 

1011011 

133 

Puncture  table 

0111110111111 

07677 

1111111111111 

17777 

5.  SIMULATION 

The  capabilities  of  a  RCPC  code  in  a  binary  symmetric 
channel  (BSC)  with  additive  white  Gaussian  noise  (AWGN)  was 
determined  through  a  Monte  Carlo  computer  simulation.  This 
result  was  then  compared  with  the  calculated  upper  bound.  The 
communications  system  model  for  this  simulation  is  given  in 
Fig.  2.  The  model  is  divided  into  three  major  sections, 
transmitter,  communication  link,  and  receiver.  These  section 
are  described  in  the  following  paragraphs. 
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TRANSMITTER 


Fig.  2.  Simplified  block  diagram  of  communication 

model . 
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Transmitter 


The  transmitter  is  modelled  as  a  random  binary 
source,  a  convolutional  encoder,  and  code  puncturer.  Random 
digit  data  is  obtained  by  rounding  to  the  nearest  binary 
integer  the  output  of  a  random  number  generator.  The  number 
generator  is  of  the  multiplicative  congruential  type  which 
generates  numbers  uniformly  distributed  between  0  and  1.  The 
encoder  is  a  rate  one-half  encoder  with  a  constraint  length  and 
generator  determined  by  the  user.  The  code  puncturer  is  as 
shown  in  Fig.  1,  with  the  period  and  puncturing  table 
determined  by  the  user. 

Communication  Link 

The  communication  link  is  modelled  as  a  binary  bit 
stream  with  AWGN  causing  individual  bits  to  be  corrupted.  The 
determination  of  whether  to  corrupt  a  bit  is  determined  by  the 
theoretical  probability  of  bit  error  due  to  noise  (Pb[noise]), 
which  is  an  input  to  the  simulation.  A  separate  random  number 
is  created  by  the  random  number  generator  described  above.  If 
the  random  number  is  less  than  or  equal  to  Pb[noise],  the  value 
of  the  bit  is  inverted  (0  changed  to  1  or  1  changed  to  0).  If 
the  random  number  is  greater  than  P^l noise J,  the  value  of  the 
bit  is  unchanged. 

Receiver 


The  receiver  is  modelled  as  a  hard-decision 
maximum-likelihood  Viterbi  decoder  that  has  been  modified  to 
decode  RCPC  codes.  The  decoded  data  is  stored  for  later 
analysis . 
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6 .  RESULTS 


The  simulation  of  the  communication  system  model  was 

used  to  determine  the  accuracy  of  the  bound  on  the  probability 

of  bit  error.  Using  the  same  generator  matrix  and  puncturing 

table  as  were  calculated  in  section  3,  simulations  were  run  to 

_2 

generate  output  bit  error  rates  (BER)  for  input  P^  of  10  to 
10”®.  In  addition,  simulations  were  run  for  an  unpunctured  code 
of  the  same  generator  matrix.  These  results  are  shown  in  Fig. 

3.  For  the  range  of  interest,  and  specifically  a  BER  of  10”®, 
the  upper  bound  and  the  simulation  results  are  relatively 
close . 

The  results  shown  in  Fig.  3  also  demonstrates  the 
performance  degradation  caused  by  using  the  punctured  code 
instead  of  the  unpunctured  code.  The  input  level  required  to 
produce  a  BER  of  10”5  are  given  in  Table  2. 


Code 

Type 

Input  P^  for  an  c 

}utput  BER  •  10  ® 

Upper  Bound 

Simulation 

RCPC  code 

0.0133 

0.0148 

Convolutional 

code 

0.0171 

0.0192 

Table  2.  Input  P^  required  for  an  output  BER  of  10”® 
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Viterbi  decoding 

-  Punctured  bound,  rate  =  13/24 

X  Simulation,  rate  =  13/24 


Comparison  of  simulation  and  bounds 


To  quantify  the  comparison  of  RCPC  codes  to  standard 

convolution  codes  the  modulation  scheme  must  be  taken  into 

consideration.  The  modulation  scheme  used  in  this  paper  is 

M-ary  frequency  shift  keying  (MFSK),  where  M-8.  From  Proakis 

[4],  the  probability  of  symbol  error  (P  )  is  given  by 

s 


M-l 

M_l'|  -nkrb/(n+l) 

.  n  J  n+1  * 

n-1 


...  (1) 


where  rb  -  signal  to  noise  ratio  (SNR)  per  bit 
k  ■  3,  number  of  bits  per  symbol. 


The  bit  error  probability  (P^)  is  obtained  from  Pg  by 

„k-l 


2k  -  1 


..  (2) 


By  combining  equations  (1)  and  (2)  we  get  Pfa  expressed  as  a 
function  of  yb,  designated  as 


Pb  -  /[rbJ  ...  (3) 

The  value  of  rb»  given  Pb,  is  determined  using  Newton's 
iteration  method  from  equation  (3). 


The  upper  bounds  for  the  RCPC  code  and  the  equivalent 
convolution  code  are  shown  in  Fig.  4  versus  the  signal  to  noise 
ratio  for  MFSK.  For  reference,  the  uncoded  case  is  also  shown 
in  Fig.  4.  The  degradation  in  using  the  RCPC  code  over  the 
normal  convolution  code  is  0.3  dB.  The  advantage  gained  from 
using  a  RCPC  code  is  that  the  link  is  now  able  to  handle  a 
variable  rate. 
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_  Unpunctured,  code  rate**l/2 


Comparison  of  error  bounds. 


7.  IMPROVEMENT  OF  TIMING  ROBUSTNESS 

The  solution  presented  in  the  previous  sections 
treats  the  control  bits  as  requiring  the  same  protection 
against  jamming  and  noise  as  the  data.  An  error  in  a  data  bit 
would  not  propagate  to  other  bits,  but  an  error  in  a  control 
bit  would  cause  the  timing  of  all  ensuing  bits  to  slide  by  at 
least  one  bit.  In  an  military  operational  environment  this 
would  be  a  serious  error,  causing  the  decryption  and  decoding 
to  lose  synchronisation.  The  link  would  have  to  be 
re-synchronised. 

A  method  to  increase  the  resilience  of  the  control 
bit  to  noise  and/or  jamming  is  to  add  redundant  bits.  A  scheme 
using  a  27/48  RCPC  is  proposed  in  Table  2.  This  scheme 
increases  the  Hamming  distance  (a  measure  of  the  error 
correction  capability)  between  critical  control  words,  thereby 
decreasing  the  probability  of  error.  Bit  positions  1  to  4 
contain  the  control  bits  and  indicate  whether  there  are  23  or 
25  bits  of  data.  If  there  are  25  bits  of  data,  bit  positions  1 
to  4  also  determine  the  value  of  data  bits  24  and  25.  Bit 
positions  5  to  27  contain  data  bits  1  to  23. 


Table  2.  Definition  of  control  bits  for  jammer  resistance. 
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The  two  critical  control  words  in  the  proposed  scheme 
are  the  ones  which  determine  if  there  is  23  bits  of  data  or  25 
bits  of  data.  The  Hamming  distance  between  these  two  control 
words  is  3,  indicating  that  the  scheme  can  correct  single  bit 
errors.  Therefore,  for  an  timing  error  to  occur,  two  of  the 
control  bits  must  be  in  error.  This  improvement  in  timing 
robustness  is  demonstrated  in  the  comparisons  in  table  3. 


Probability 
of  data 
error 

l.OxlO-3 

l.OxlO-4 

l.OxlO-5 

l.OxlO-6 

Probability 
of  timing 
error 

6 . OxlO-6 

6. OxlO-8 

6. OxlO-10 

6. OxlO-12 

Table  3.  Comparison  of  probability  of  error  for  data  and 
timing 

To  implement  the  scheme  defined  above  the  number  of 
data  bits  plus  control  bits  must  be  increased  to  27.  In  order 
for  the  overall  coded  data  rate  to  remain  constant,  a  rate 
27/48  (or  9/16)  code  will  be  required.  The  performance  of  this 
code  is  compared  with  the  previous  codes  in  Fig.  5,  showing  a 
degradation  of  0.5  dB  from  the  unpunctured  code.  If  the 
allowable  variation  in  the  data  rate  was  reduced  to  ±2%,  a 
scheme  could  be  implemented  with  a  rate  51/96  (or  17/32)  code. 
The  performance  of  this  code  is  also  presented  in  Fig.  5, 
showing  a  degradation  of  only  0.2  dB  from  the  unpunctured  code. 

8.  OTHER  USES  OF  RCPC 

The  previous  sections  have  described  how  a  RCPC  can 
be  used  to  implement  a  concept  to  support  equipment  with 
non-synchronised  clocks.  This  section  will  describe  how  a  RCPC 
can  be  used  to  transmit  data  rates  which  are  not  exactly 
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Upper  bounds 

-  Unpunctured,  code  rate=l/2 


Comparison  of  error  bounds 


matched  to  the  data  rate  of  the  data  communication  system, 
without  using  filler  bits. 


Many  MILSATCOM  systems  assign  communication  resources 
in  increments  of  1200  bps.  Therefore,  a  16  kbps  digital  voice 
transmission  would  normally  require  a  16.8  kbps  link,  using 
filler  bits.  Using  a  rate  20/39  code,  the  16  kbps  voice  could 
be  transmitted  over  a  nominal  15.6  kbps  link,  requiring  no 
filler  bits  and  using  less  bandwidth. 

A  Tl  carrier  is  defined  as  1,536  kbps  of  data  and 
8  kbps  of  framing,  for  an  overall  data  rate  of  1,544  kbps.  A 
typical  communication  system  handling  medium  data  rate 
information  would  use  transmission  rates  of  32  kbps,  or 
multiples  of  32  kbps.  This  would  imply  that  a  1,544  kbps  data 
rate  would  require  a  1,568  kbps  channel  using  filler  bits. 
Instead  of  using  this  approach  a  rate  25/48  RCPC  would  allow 
the  full  Tl  carrier  (1,544  kbps)  to  be  transmitted  over  the 
nominal  1,536  kbps  link  without  filler  bits. 


9.  CONCLUSION 


The  concept  described  in  this  document  provides  a 
simple  method  of  implementing  a  scheme  to  allow  a  variation  of 
±4%  in  the  source  data  rate  for  a  system  with  a  fixed  symbol 
transmission  rate.  The  degradation  in  performance  when  compared 
to  the  existing  scheme,  which  allows  no  variation  in  the 
transmission  rate,  is  0.3  dB.  This  is  considered  a  very 
acceptable  performance  sacrifice.  In  addition,  little 
modification  is  needed  to  the  data  link  standard  and  technical 
complications  are  kept  to  a  minimum. 
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APPENDIX  A 


Calculation  of  the  Upper  Bound 
For  Convolutional  Codes 


This  appendix  describes  a  method  to  determine  the 
performance  bounds  on  specific  convolutional  codes  using  binary 
symmetric  channels.  The  initial  work  was  given  in  [3]  for 
standard  convolutional  codes  with  extensions  in  [2]  for  RCPC 
codes . 


The  decoding  of  a  convolutional  code  is  usually 
described  as  moving  through  a  Trellis  diagram,  as  shown  in 
Fig.  Al.  If  the  decoder  employs  the  Viterbi  algorithm  it  is 
classed  as  maximum-likelihood  decoder  in  that  it  selects  the 
path  with  the  maximum  likelihood  function  (or  metric)  as  the 
correct  path  in  the  Trellis.  As  shown  in  Fig.  Al  the  path 
chosen  is  the  path  with  the  largest  metric.  The  measure  of  the 
performance  of  a  code  is  usually  given  by  the  bit  error 
probability,  P^.  In  convolutional  codes,  and  RCPCs,  it  is  not 
possible  to  explicitly  define  the  probability  of  bit  error 
(P^),  but  P^  can  be  upper  bounded. 

Before  calculating  the  upper  bound  on  P^,  we  will 
turn  to  the  determination  of  the  more  readily  available  P  ,  the 
error  probability  per  node.  Referring  to  Fig.  A2,  we  see  two 
paths  originating  at  node  j.  Without  loss  of  generality,  we 
take  the  all  zeros  path  to  be  the  correct  path  and  the  other 
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Fig.  A2  .  Example  of  an  incorectly  decoded  path  in  a  Trellis  diagram. 


path  to  be  an  incorrect  one  chosen  by  the  decoder.  For  this 
error  to  occur  the  correct  path  over  the  unmerged  segments  must 
have  a  lower  metric  than  the  incorrect  path.  Therefore,  we  can 
upper  bound  P  (employing  the  union  bound)  by 


Pn(  j )  <  ^  Pr  [  AM(Xj',Xj)  >  0  ) 

x j '  is  an  element  of  X'(j) 

where  x ^ *  -  is  the  incorrect  path  diverging  at, 
node  j 

X'(j)  -  is  the  set  of  all  such  paths,  and 


. . . (Al) 


AM( x j ' , Xj )  -  is  the  metric  difference  between 
the  correct  and  incorrect  path. 

Each  term  of  the  summation  is  the  pairwise  error  probability 
for  the  two  code  vectors  over  the  unmerged  segment  which  yields 


Pn(j)  S  Y.  a(d)  Pd 


where  d^  *  minimum  free  distance  of  code, 

a(d)  ■  incorrect  path's  Hamming  distance  d  from 

the  correct  path,  and 
Pj  -  pairwise  error  probability  for  code 
vectors  at  distance  d. 


(A2) 


The  value  P^  can  now  be  bounded  by  weighting  each 
term  of  the  union  bound  by  the  number  of  bit  errors  which  occur 
on  that  path.  This  corresponds  to  the  number  of  "l"s  in  the 
data  sequence  over  the  unmerged  segment  of  the  incorrect  path 
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(the  all  zeros  path  was  selected  as  the  correct  path). 
Therefore , 


pb  i  Y.  E  ia<d'i)E,d 

i-l  d-df 

where  i  *  number  of  "l"s  in  the  data  segment  of 
the  incorrect  path. 


...  (A3) 


The  exact  values  of  P^  can  be  determined  for  a  binary  symmetric 
channel  (BSC)  using 


X  (k)  pk  (1"P)d“k  d  odd 

k-(d+l)/2 

“d/2  Pd/2d-P)d/2  +X  (k)  rk,1-P>d"k  d  even  •••  <A4> 

k-(d+l)/2 
where  p  -  input  error  probability. 

To  use  the  above  method  for  RCPC  codes  the  time 
varying  nature  of  the  code  must  be  accounted  for.  P  (period  of 
punctuation)  different  starting  points  for  diverging  paths  need 
to  be  considered.  Equation  (A3)  then  becomes 
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.  ..  (A5) 


pb  i  ?  E  E  ia(d'i)pa 

i-l  d-df 

The  values  a(d,i)  and  i  were  determined  from  a 

computer  program  that  searched  for  all  possible  paths  up  to 

length  20.  Path  lengths  greater  than  20  had  very  little  effect 

2 

on  the  value  of  the  upper  bound  for  input  probabilities  <  10-  . 
This  program  works  with  all  rate  1/2  convolutional  codes  up  to 
constraint  length  8  or  any  RCPC  derived  from  these 
convolutional  codes.  Using  a  convolutional  code  of  rate  1/2, 
constraint  length  7,  and  generator  matrix  171g,  133Q,  the  upper 
bound  is 

pb  -  36Pd(10)  +  21lPd(12)  +  1404Pd(14)  +  11633Pd(16) 

+  77433Pd(18)  +  502690Pd(20) . 

This  is  the  same  result  obtained  by  Marvin  K.  Simon  et  al  [5]. 
The  upper  bound  for  the  optimised  RCPC  code  using  the  above 
convolutional  code  and  puncturing  values  of  07677g  and  177777g 
is 

pb  "  3Pb(8)  +  26pb(9)  +  63pb(10)  +  162Pb(ll)  +  457Pfa(12) 

+  1298Pb(13)  +  3908Pb(14)  +  9896Pb(15)  +  22465Pb(16) 

+  48068Pfa(17)  +  95527Pb(18)  +  177398Pb(19)  +  302119Pb( 20 ) . 
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